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Abstract. The phase space corresponding to a particular four-parameter family of initial data for 
the gravitational collapse of a spherically symmetric dust cloud is investigated. In a certain Umit 
of the parameters, this family reproduces the case of homogenous mass density -constant mass 
distribution- and zero initial velocity, while in another limit, it generates a globally naked singularity. 
We show that for initial data characterizing flat density profiles, as well as large initial velocities, 
the probabihty of fomning a globally naked singularity is low. 
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INTRODUCTION 

In General Relativity, the singularity theorems [1] predict that a spacetime singularity 
forms when a sufficiently large mass is concentrated in a small region of an asympto- 
tically flat spacetime, as for instance occurs in the complete gravitational collapse of a 
star. An open question is whether this singularity is naked (light rays emanating from 
it are visible by distant observers), or hidden (a black hole censures the singularity so 
that no information escaping from it can be detected outside of the event horizon). Ac- 
cording to the weak cosmic censorship conjecture [2], the second must occur in the 
generic, realistic matter case [3]. In this work, we consider the simple model of a spheri- 
cally symmetric dust cloud, whose collapse can lead to the formation of globally visible 
shell-focusing singularities [4, 5, 6, 7] without fine-tuning the initial data, indicating that 
they are stable within this particular model. Although the model is restricted to spherical 
symmetry and a pressure-less fluid (and thus its matter content is not fully realistic), 
the formation of naked singularities is interesting by itself. In this article, we construct 
a four-parameter family of initial data for the spherical dust collapse. For certain val- 
ues of the parameters, it gives a hidden singularity from constant mass density and zero 
initial velocity, while for other values, it gives a naked singularity. We construct the cor- 
responding conformal diagrams using numerical techniques [6, 8] and analyze the phase 
space of the parameters of the initial data. 

TOLMAN-BONDI DUST COLLAPSE MODEL 

Inside a self-gravitating, spherically symmetric collapsing dust cloud of finite radius R\, 
in co-moving, synchronous coordinates {x,R) [9], the metric g and four- velocity u are 



g=-dz^+ dR^ + r{z,R)\di»^ + sm^i»d(p^), u = |^ , (1) 

where, for given T, r(T,i?) is the areal radius of a matter shell of initial radius R, this 
means that the coordinate R is such that r(0,i?) = R. Our assumptions below exclude 
shell crossing singularities; therefore, each matter shell is uniquely characterized by R. 
Here and in the following, a prime denotes partial derivation with respect to R and a dot 
with respect to T. From Einstein's equations, the time evolution of r is determined by 

^-riT,Rf + ViriT,R),R)=EiR), y(r,i?) := (2) 

where m{R) is the Misner-Sharp mass function [10] and E{R) = vo{R)^/2 — m{R) /R is 
the total energy for each shell R, given in terms of the initial data, which consist of the 
profiles for the initial velocity vo(^) and the initial density Po{R). Now, it is convenient 
to introduce the new functions 



Note that c{R) is directly proportional to the mean density and q{R)^ is the ratio between 
the total and initial potential energy. The initial data must satisfy specific conditions such 
as regularity, smoothness, etc. (see Ref. [6] for details and explanations), which imply 

c{R)>0, c'{R)<0, 0<q{R)<l, q'{R)>0, (3) 

as well as q'{R)/R > whenever c'{R)/R = 0, for R>Q. This guaranties [6] the 
existence of light rays escaping from the central singularity, making it visible, at least 
to local observers. These assumptions do not cover the Oppenheimer-Snyder case [11], 
(£■ = and constant mass density profile), however, it can be recovered by approximation 
from data satisfying our conditions. 

INITIAL DATA AND CONFORMAL DIAGRAIMS 

A particular choice for the functions c and q is the four-parameter family 

c(7?):=co 1-^^f-j , ?W:=?o + ?if^J , 0<i?<i?i, (4) 

where /?i > is the initial radius of the cloud (for R > Ri we have Schwarzschild 

spacetime). The parameters are subject to n > 1, < R^cq < (3 + 2n) /In, < qo < 
1, and < qi < I — qo, which guarantee the satisfaction of conditions (3) on the 
interval [0,7?i]. The corresponding initial density and velocity profiles are (po,vo) = 
{{R^cy/{SnGR^), -Ry/ (1 -q^)c). Explicitly, the resulting density profile is 



Po(R) = 



1 



R 



0<R<Ri. 



The important feature of this choice is that po becomes constant as n — )■ oo, moreover, the 
initial velocity becomes zero as qo \ and qi ^ (Note that the Oppenheimer- Snyder 
case is recovered when both qo,qi 0). Figures 1 and 2 (left panel) are examples of 
numerically generated conformal diagrams - in coordinates {T,X), in which the radial 
part of the metric (1) takes the form Q.{T,X)^ {—dT^ + dX^^, with a conformal factor 
^2(r,X) > -, with initial data given by (4) with n = Al (nearly constant density), 
qQ — 0.98 and q\ = 0.01 (nearly zero initial velocity) for both diagrams. The only 
parameter that has a different value in Figs. 1 and 2 is cq. In Fig. 1, cq = 0.063, so 
the central density po = 3cq/{%tzG) is sufficiently small to generate a globally naked 
singularity, as predicted by theorem 2 in Ref. [6]. In Fig. 2, cq has the limit value 
Co = 1.0, which gives a hidden singularity. Its right panel shows the exact diagram 
for homogeneous -constant- density and zero initial velocity. The differences between 
the two panels are due to the implementation of different boundary conditions for their 
construction. Note that, while Fig. 1 shows the complete spacetime, including the outside 
of the collapsing cloud, for simplicity Fig. 2 only shows the interior spacetime of the 
cloud. Also note that in Fig. 1 the apparent horizon in spacelike everywhere inside the 




FIGURE 1. Conformal diagram for the Tolman-Bondi dust collapse with initial data given by Eq. (4), 
with = 0.98, qi = 0.01, n = 41.0 and cq = 0.063. The notation "AH", "EH" and "CH" refer to the 
apparent, event and Cauchy horizons, respectively. In this case, a portion of the singularity is globally 
naked. The dotted line corresponds to the dust particle trajectory with initial areal radius Rq = 0.957?i. 
The diagram was generated following the numerical algorithm described in Ref. [6]. 

cloud, while in Fig. 2 it has a timeUke portion. 



PHASE SPACE 

The phase space corresponding to the four-parameter family of initial data given in 
Eq. (4) is studied in this section. The left panel of Fig. 3 shows a subset of this space for 
fixed q. The critical -continuous- line in the n-co-plane divides the regions corresponding 
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FIGURE 2. Left panel: Conformal diagram for the Tolman-Bondi dust collapse with initial data given 
by Eq. (4), with the same parameter choice as in Fig. 1, except for cq = 1.0. This is a hidden singularity 
inside a black hole. This singularity has a null portion as in Fig. 1, but here it is too small to be seen without 
a zoom. This is an approximation to the scenario of uniform density and zero initial velocity. Right panel: 
Exact conformal diagram for the case of uniform density and zero initial velocity, see Ref . [6] . 



to initial data giving rise to black holes and globally naked singularities (shaded region). 
It results that flat density profiles (large n) favor the formation of black holes. It is also 
inferred that the critical value for cq becomes asymptotically constant for large n. On 
the other hand, the dashed line corresponds to the upper bound for cq obtained from 
theorem 2 in Ref. [6]. Note that although this bound is not optimal, it describes the 
correct qualitative behavior for the critical line. The right panel of Fig. 3 shows the 
critical line in the n-CQ-plane for different values of qo and qi = 0.01. Since 1 — is 
the ratio between the initial kinetic and potential energies, this figure shows that large 
initial velocities in the negative radial direction favor the formation of black holes. All of 
this is consistent with the conclusions obtained in section 5.2 of Ref. [6] using a different 
family of initial data. Nevertheless, an important new feature obtained from Fig. 3 is that, 
for large initial velocities, the region for naked singularities becomes small for n >2, 
which means that it becomes very difficult to obtain a naked singularity in this case. 

CONCLUSIONS 

The Tolman-Bondi model for the gravitational collapse of a spherically symmetric dust 
cloud with fairly general initial density and velocity profiles was considered. A four- 
parameter family of initial data was constructed such that in a certain limit of the 
parameters, it arbitrarily approaches the case of homogenous -constant- density with 
zero initial velocity, while in another limit, a globally naked singularity is obtained. 
Some important remarks are: (i) With an arbitrarily flat density profile, it is always 
possible to obtain a globally naked singularity for small central density, (ii) Just as in 
the case of constant density and zero initial velocity, a portion of the apparent horizon 
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FIGURE 3. Left panel: A cut of the parameter space of the four-parameter family of initial data (4), 
corresponding to the subset with fixed qo = 0.98 and qi = 0.01. The shaded region represents initial data 
giving rise to globally visible singularities. The dashed line describes the upper bound for cq given from 
theorem 2 in Ref. [6]. Right panel: The critical line, dividing the regions of initial data giving rise to black 
holes and globally naked singularities, respectively, for fix qi — 0.01, and different values of ^o- Note that 
the region for naked singularities becomes very small for large initial velocities and flat density profiles. 

may be timelike. (///) Flat density profiles favor the formation of black holes, (iv) Large 
initial velocities favor the formation of black holes, (v) For flat density profiles and large 
initial velocities, the probability of obtaining a naked singularity is very low. Since at the 
late stage of a collapse the velocity is high, this means that for a sufficiently flat profile 
it is difficult to obtain a naked singularity. 
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